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We study the effects of instantons on the charged pion electro-magnetic formfactor at intermedi- 
ate momenta. In the Single Instanton Approximation (SI A), we predict the pion formfactor in the 
kinematic region Q 2 — 2 — 15 GeV 2 . By developing the calculation in a mixed time-momentum rep- 
resentation, it is possible to maximally reduce the model dependence and to calculate the formfactor 
directly. We find the intriguing result that the SIA calculation coincides with the vector dominance 
monopole form, up to surprisingly high momentum transfer Q 2 ~ 10 GeV 2 . This suggests that 
vector dominance for the pion holds beyond low energy nuclear physics. 
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I. INTRODUCTION 

Bridging the gap between the non-perturbative and the 
perturbative sector of QCD is a central step toward our 
understanding of the strong interaction. In this context, 
the electro-magnetic formfactor of the charged pion 
is of great interest. It is, at low momenta, extremely well 
reproduced by the vector dominance model. In addition, 
at very high momentum transfer, it is constrained by per- 
turbative (p)QCD predictions (for a review on hadronic 
formfactors, see e.g. [Q). The asymptotic behavior for 
large space- like momentum transfer, Q 2 = — (p—p') 2 > 0, 
is derived in a closed form in perturbation theory [g, |^, |J , 



Q 2 F V (Q 2 ) 



16 tt f 2 a s (Q), 



(1) 



where /„■ = 92.4 MeV denotes the pion decay constant. 

A comparison of the asymptotic behavior and the ex- 
perimental data determines the momentum scale where 
the perturbative regime of QCD is reached. Recently, 
the charged pion formfactor has been measured very 
accurately at momentum transfers 0.6 GeV 2 < Q 2 < 
1.6 GeV 2 by the Jefferson Laboratory (JLAB) F n col- 
laboration j5j and lead to quite surprising results. Not 
only are the data at highest experimentally accessible 
momenta still very far from the asymptotic limit, but 
the trend is away from the pQCD prediction (see Fig. [I]) . 

Moreover, it is quite remarkable that the data are 
still completely consistent with the vector dominance 
monopole fit, 
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Mp = 770 MeV, 



(2) 



at relatively high momentum transfer (Q 2 w lGeV 2 ). 
Clearly, the measurements currently undertaken at JLAB 
in the region 0.5 GeV 2 < Q 2 < 6 GeV 2 are very much 
needed and will provide information about whether the 
perturbative region is reached at that scale. 

The charged pion formfactor has attracted a lot of at- 
tention from the theoretical side @, § §, [k| 0, [jj, H)- 
Despite this debate, we feel that there are still a number 
of open question. What are the leading non-perturbative 
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FIG. 1: The recent JLAB data for Q Fn{Q ) in comparison 
with the asymptotic pQCD prediction (thick bar, for a typical 
a s ~ 0.2 — 0.4 in Eq. (p])), the monopole fit (dashed line), and 
our SIA calculation (solid line). The SIA calculation is not 
reliable below Q 2 ~ lGeV 2 . The solid circles denote the 
SLAC data 101. 



effects responsible for the deviation from pQCD, Eq. (jl]), 
at intermediate momenta? Where can we expect the 
transition to pQCD, is it within experimentally accessible 
momentum transfers? Is there a microscopic explanation 
of the success of the monopole form? Can it be justified 
at so large momenta, where the vector dominance model 
should be inadequate? 

In this work we suggest answers to these questions. In 
particular, we study the effects of the leading (i.e. zero- 
mode) interaction of light quarks with the intense classi- 
cal vacuum fields (instantons) on the pion formfactor, for 
Q 2 > lGeV 2 . We show that, in this kinematic region, 
the pion formfactor is dominated by the interaction of 
quarks with a single instanton. 

The pion plays a special role among the hadrons. Not 
only is it nearly massless (which is explained by the Gold- 
stone theorem), but unusually compact as well. Phe- 
nomenologically, this is seen e.g., from the rather large 
electro-magnetic mass splittings between the charged and 
neutral pion. Theoretically, it was explained by studying 
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instanton-induced effects on the two-point correlators, 
see |l4| for a review. In the pion channel (as well as other 
scalar and pseudoscalar channels) the instanton contribu- 
tion to the two-point correlator can be represented by the 
zero-mode terms in the quark propagators, and as a re- 
sult the effect is enhanced (relative to e.g., the vector or 
axial channels) by a factor l/(m*p) 2 , where p denotes the 
typical instanton size, and m* is an effective parameter 
with the dimension of a mass, defined and discussed in 
detail in . Numerically, such an enhancement factor is 
about 30, and paramctrically it is the inverse dilutcncss 
of the instanton ensemble. Due to the presence of such a 
large factor, instanton-induced forces become dominant 
in the pion pseudoscalar correlator, starting from rather 
small distances. The same factor is present in the in- 
stanton contribution to the three-point function, which 
is related to the pion formfactor. Although, the final re- 
sults we obtain are independent of the value of m*, the 
window of applicability of our method does depend on it. 

This feature, however, is not generically related to 
the pion itself and depends on the particular three-point 
function under investigation. For example, the enhance- 
ment is absent, when one considers the pion contribution 
to the axial correlator. Similarly, there is no such en- 
hancement of the 77*7T° neutral pion transition formfac- 
tor. The relevant instanton effects for this process are not 
due to (enhanced) zero modes, but are either related to 
non-zero mode propagators in the instanton background 
or to multi- instanton effects, which are suppressed by the 
instanton dilutcncss. This conclusion is nicely supported 
by recent CLEO measurements of this formfactor, which 
indeed show that the asymptotic pQCD regime is reached 
much earlier, at Q 2 ~ 2GeV 2 [0. 

The first calculation of the pion formfactor in the SIA 
was performed by Forkel and Nielsen |Tij ], who comple- 
mented QCD sum rules by the instanton-induced contri- 
bution. In such an approach, however, a model descrip- 
tion of the continuum of excitations with the quantum 
numbers of the pion is needed, in order to connect the 
electro-magnetic three-point function to the pion form- 
factor . In order to avoid such additional model de- 
pendence, Blotz and Shuryak proposed a different ap- 
proach |19| , in which large-sized three-point functions 
obtained both from simulations in the Instanton Liq- 
uid Model (ILM) and from its spectral decomposition 
(which depends on the pion formfactor) were compared 
directly in coordinate space. This analysis revealed that, 
at large distances, the ILM results were completely con- 
sistent with the monopole fit. 

In this work we follow a different approach, which 
is based on the mixed time-momentum representation 



widely used in lattice calculations (see e.g., (20|). We 
show that this way it is possible to isolate the pion pole 
and yet make predictions in momentum space (and there- 
fore compare directly with experiments), without any 
model of the continuum contribution. 

This approach has several other advantages. Firstly, 
the model dependence is maximally reduced: Our re- 
sults depend only on one dimensional parameter, the av- 
erage instanton size, p ~ 1/3 fm, which was proposed 
many years ago JL8| and checked a number of times by 
phenomenological studies |l4| as well as lattice simula- 
tions |2l], |2^, |2^, [24|, |2f| . The calculation is independent 
of the properties of the instanton liquid ensemble such 
as the instanton-instanton interaction and the instanton 
density. Moreover, we will argue that the physical con- 
tent of the calculation is more transparent and the dy- 
namical effects of the instantons on the propagation of 
the quarks can be seen naturally. 

The central prediction following from our SIA calcula- 
tion is shown in Fig. [I] in comparison to the recent JLAB 
measurements. 

The main physical point of our analysis is that, if the 
strongly attractive 't Hooft interaction is taken into ac- 
count, the asymptotic perturbative regime of the pion 
formfactor is reached much later than in other models, 
which do not include this force. Such a qualitative dif- 
ference implies that the upcoming JLAB measurements 
will provide an important test of instanton forces in the 
pion. We also find the very intriguing result that the 
instanton contribution to the formfactor is completely 
consistent with the monopole form at intermediate mo- 
mentum transfers, 2GeV 2 < Q 2 < lOGeV 2 , where the 
vector dominance model has no justification. It is a first 
microscopic study showing that all other resonances with 
quantum numbers of the p (except the p itself) are not 
seen in the pion formfactor. For large momentum trans- 
fers, Q 2 > 20GeV 2 , our SIA breaks down, as it is nec- 
essary to increase the distances in order to isolate the 
pion ground state. At these needed distances, however, 
the correlation functions will become sensitive to multi- 
instanton effects. 



II. INSTANTON CONTRIBUTION TO THE 
FORMFACTOR 

We consider the spatial Fourier transforms of the Eu- 
clidean three-point function and two-point function, 

G„(t,p + q;-t,p) = J d 3 xd 3 ye- ip x+l(p+q) y 

x (0| i 5 (t, y) J M (0, 0) x) |0), (3) 



1 Notice that in |17| Forkel and Nielsen used the numerical value 
for the effective parameter m* , which was derived in the mean 
field approximation |fi| . However later a significantly smaller 
value (and larger enhancement) has been derived from numerical 
simulations of the Instanton Liquid Model |15| . 



G(2<,p) = |d 3 xe l P- x (0|j 5 (i,x)jt(-i,O)|0), (4) 

where the pseudo-scalar current j§(x) = u(x) 75 d(x) ex- 
cites states with the quantum numbers of the pion and 
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Jfj,(0) denotes the electro- magnetic current operator. In 
the large t limit (at fixed momenta), both correlation 
functions are dominated by the pion pole contribution 
and the ratio of the three-point function to the two- 
point function becomes proportional to the pion form- 
factor 20 . In the Breit frame, p = — q/2 and Q 2 = q 2 , 



one has simply 



G 4 (f,q/2;-t,-q/2) 
G(2t,q/2) 



F^{Q 2 ). 



(5) 



Notice that the LHS of Eq. (||) should not depend on t, 
for t large enough. Below we demonstrate that, for the 
pion, this is achieved already for t ~ 0.6 fm. 

We will evaluate the LHS of Eq. (||) in the SIA. In prin- 
ciple, it is not obvious that such an approach is justified, 
as the mean Euclidean distance between two instantons 
in the vacuum is about 1 fm. Therefore, if t ~ 1 fm, one 
would expect many instanton effects to play a non negli- 
gible role. However, two of the authors showed that the 
pion and nucleon three-point functions, evaluated in the 
SIA, agree with the results of numerical simulations in 
the instanton liquid model up to distances of the order 
of 1 fm Moreover, they found that the ratio of three- 
to two-point function is dominated by single instanton ef- 
fects to even larger distances [^6| . This result enables us 
to reliably evaluate this ratio in the SIA up to Euclidean 
distances of ~ 1.4 fm. A possible physical interpretation 
of this fact is that in order for the scattered parton to 
recombine in the same final state (and therefore have an 
elastic formfactor), it is sufficient that the parton scatters 
off a single instanton during the process. 

We expect that the main contribution to the time- 
momentum correlation functions comes from distances of 
the order of the inverse conjugate momenta, |x| < l/|p| 
an d |y| ;$ l/|p + q|. Therefore, in the Breit frame and 
with the choice t ~ 0.7 fm, one obtains that the corre- 
lation functions under consideration can be evaluated in 
the SIA, for momenta |q| > 0.6 — 1 GeV. We note that 
one of the authors applied the same method to evaluate 
the pion and nucleon dispersion curves f27j| . Agreement 
between the SIA calculation and experiment was found 
for t > 0.7 fm and |p| > 1 GeV. 

Having assessed the applicability of our approximation, 
we proceed to the calculation. The connected three-point 
and two-point functions are given by 2 

G 4 (t, p + q; -t, p) = (e„ - e d ) 

x j d 3 xd 3 ye _,;p ' x+,;(p+q) ' y (Tr{7 5 5(t,y;0,O) 

x 74 S(0, 0; -t, x) 75 S(-t, x; t, y)}), (6) 



2 A discussion of the contribution from disconnected diagrams can 
be found in jis[ ]. 



G(2t,p) = J d 3 xe lp x (Tr{ 75 S*(t,x;-<,0) 

x 75 S(-t,0;i,x)}), (7) 

where S(y, x) denotes the quark propagator in the back- 
ground of a gauge field, the Tr is over spin and color, and 
the brackets ( • } denote the average over all gauge field 
configurations. 

We express Eqs. (|^) and (j^) in terms of "wall-to-wall" 
(W2W) correlators, defined as the spatial Fourier trans- 
forms of the quark propagators 

S(t',p';t,p) = /'d 3 xd 3 ye ip '- y - ip - x S'(y,x). (8) 

This is achieved by insertions of appropriate delta func- 
tions at each vertex 

f d 3 k f d 3 l f d 3 xn 

x J ^(Tr{ 7 g5(t,k;0 ) m + n)74S'(0,m;-t,l) 

x 75 S'(-t,l-p;t,k-p-q)}) ) (9) 



G(2t,p) 



d 3 k 



d 3 l 



(2tt) 3 J (2tt) 3 J (2?r) 



d 3 ! 



x (Tr{7 5 5(t,k;-t,l) 75 5(-t,m;t,k-p)}). (10) 

These W2W correlation functions are completely gen- 
eral. Their graphical representation, as an example for 
the three-point function, Eq. (||), is given in Fig. || (A). 

A model dependence is introduced, when we evaluate 
the gauge field average in the SIA. This is achieved by 
taking the quark propagators in the background of an 
instanton and averaging over the collective coordinates 
(instanton position, size and color orientation). The in- 
tegration over the color orientation is trivial, while the 
integration over the instanton size is weighted by a distri- 
bution function. We use a distribution for the instanton 
size which takes into account the small-size limit. The 
latter has been calculated by 't Hooft, by considering 
perturbative fluctuations around the instanton solution 
For large-sized instantons, however, the distri- 
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bution has to be cut off. We use a Fermi distribution 
function to account for the suppression of large-sized in- 
stantons, where we fit the width and the mean instanton 
size to lattice QCD results (for a summary on lattice 
QCD results, see e.g., Negele @). This leaves for the 
single instanton density 

n{p) =nd< t Hooft (p) 77 (H) 

exp((p-p)/cr) +1 

where n denotes the average instanton density and the 
QCD small-size limit is given by 't Hooft 



't Hooft 



(P) 



JVf-5 



2N C 



xexp{-(3 2 (p) + (2N c -—) 



2b'2b[3 1 {p) 



m(/?i(p))}- (12) 
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(B) 




(D) 



FIG. 2: (A): Graphical representation of the W2W three point 
function, Eq. (^). The double lined "walls" correspond to the 
spatial Fourier integration. (B): A particular diagram con- 
tributing to (A), in lowest order pQCD. (C) and (D): The 
leading non-perturbative contributions to (A) in the SIA. The 
dashed ellipse denotes the quark (zero-mode) 't Hooft inter- 
action. Note that the momentum integrations in Eqs. ([]) 
and (|nj) ensure locality at the walls and the gauge invariance 
of the calculation. 



Here, 0i(p) and ffa{p) denote the one and two-loop beta 
functions 3 , 



Hp) 



-b\n(p Apv) 
b' 



/3 1 (p) + ^ln(21n(l/pA Py )), 
b =Y Nc ~l Nf > 

3 C 3 J Nr 



(13) 
(14) 

(15) 

(16) 



For our analysis, we use a parameterization, which sum- 
marizes recent lattice QCD results Q, p = 0.37 fm 
and a = 0.15 fm. We will also contrast this size den- 
sity to a much simpler model proposed by Shuryak |L8| , 
n(p) = n S(p — p), where p = 1/3 fm. 

The quark propagator in the field of an instanton was 
evaluated exactly in singular gauge |5lJ and consists of a 
zero-mode and a non zero- mode part, 



S ! {x,y) = S I zm (x,y) + Sj lzm (x,y). 



(17) 



The expression for S^ zm is quite involved. However, at 
small distances, \x — y\ <C p, and if the instanton is far 



away, \x — z\,\y — z\^> p 1 the non zero-mode propagator 
reduces to the free one. For this reason, most calculations 
are carried out in the zero- mode approximation (ZMA) 4 , 
in which one approximates the non zero-mode part with 
a free massless quark propagator, 



S\x, y) w S (x, y) + Sl m {x, y). 



(18) 



Let us discuss the accuracy of the ZMA for our calcu- 
lation. We recall that, in the momentum range of inter- 
est, the contributing point-to-point correlators in Eqs. (^) 
and (g), are ~ 1 fm long. It is therefore sufficient to ac- 
cess the accuracy of the ZMA for our point-to-point cor- 
relators at those distances. In Eq], two of the authors 
have evaluated several three- and two-point functions in 
the SIA with ZMA and compared them with the results 
of numerical simulations in the instanton liquid model 
(where non-zero mode effects are taken into account). 
Very good agreement was observed at distances relevant 
for our calculation. This result could either imply that 
multi-instanton and non-zero mode effects are individu- 
ally small at that scale, or that both effects are not quite 
small, but that they happen to cancel. 

In this approximation, the non-perturbative contribu- 
tions to Eq. are represented in Figs. |^ (C) and (D) and 
are straightforward, once the free and zero-mode W2W 
quark propagators are evaluated. 

The massless free W2W quark propagator is given by 



S (t',p';t,p) = (27r) s 6^(p'-p) 



,-|p| \t'-t\ 



'V)*»> ( 19 ) 



where 114 = 1 and ui — ipi/\p\, for I — 1,2,3. We note 
that the free W2W propagator is proportional to a delta 
function due to momentum conservation. The free con- 
tribution to G4 and G are then given by 



G 



(free) 



(t,P 



q;-*,pj = -5 



3 f d 3 k / k • (k - q) 



(27T)3 V |k| |k - q 
k • (k - p - q) (k - q) • (k - p - q) 



Ikllk- 



x e 



k q| |k p q| 

t(|k| + |k-q|+2|k-p-q|) ; (20) 



G (-) (2t ,p)=3/ " 



d 3 k ( k • (k - p) 



(27T) 3 V |k||k-p| 



- 1 



x e 



-2t(|k| + |k-p|)_ (21) 



The free W2W two-point function can be further simpli- 
fied, 



G (froc) (2t,p) 



3 l + 2t|p| 2t | p , 
(2tt) 2 (2i) 3 



(22) 



3 This result was carried out in the Pauli-Villars regularization 4 We note that, at large distances (\x — y\ p), the zero-mode ap- 
scheme and Apy is the corresponding scale parameter. We use proximation gives the same correlation functions as the 't Hooft 

Apy = 250 MeV. Lagrangian. 
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Next, we calculate the zero-mode W2W quark prop- 
agator. Since we are eventually interested in a gauge 
invariant quantity, and we are working in the SIA, we 
can use the instanton solution in regular gauge 5 , where 
we find the simple result: 



S'W (f , p'; t, p) = ^ f(t', p'; t, p) W 1 ^ 



(23) 



' n'- + n\ = e l (p'~p) z 



/(*',p';t,p) 



K (]v'\^{t> -z A ) 2 +p 2 ) 



x K Q (\p\ y/(t-z 4 )*+fP), (24) 

(25) 



fl 'VI 



where — (z, z&) denotes the instanton position, to* 



is 



the effective 
defined in 



arameter encoding many-instanton effects 



and 



r± - 



are color matrices. We 



remark that the spatial instanton position appears only 
coupled to the quark momentum transfer in the back- 
ground of the instanton field. This is intuitive: if the 
instanton is far away from the quark, there should be 
little momentum transfer. Due to this fact, the integra- 
tion over the spatial instanton position leads to an overall 
momentum conserving delta function. As expected, the 
momentum conservation in the SIA is a consequence of 
the translational invariance of the spatial instanton posi- 
tion. Moreover, the modified Bessel functions exponen- 
tially suppress momentum transfers much larger than the 
inverse size of the instanton, 1/p. This implies that large 
momentum transfers are driven by small instantons and 
directly probe the small-size limit of the 't Hooft measure 
for the tunneling amplitude. This argument holds for all 
W2W correlation functions evaluated in the SIA and was 
applied by one of the authors to the pion and nucleon cor- 
relators [27|. It was shown that at sufficiently large mo- 
menta, a 1 the model-dependence in such a semi-classical 
approach is removed, and the physical observables can be 
expressed in terms of Aqcd- We shall similarly study at 
what momentum transfer the pion formfactor becomes 
sensitive to the small-sized instantons only. 

The explicit evaluations of the zero-mode contribu- 
tions, Figs. H (C) and (D), lead to 



G 



(C) 
4 zm 



OC fOO 



dp 



(Iza 



d 3 k f d 3 l 



n(p) 



(2tt) 3 J (2tt) 3 

^ e-W K (\k - q| - k| £(z 4 - t)) 

x K (\ - 1 - q/2| £(z 4 + t)) K (\l\ + t))> (26) 



G 



(D) 



oo poo 



x ^ *o(|k| £(z 4 - t)) tf (|q/2 - k| £(z 4 - t)) e-W « 
x K (\ - 1 - q/2| £(z 4 + t)) Ml 1 + q| SW), (27) 



where £(x) = \J x 2 + p 2 . We have labeled the mo- 
menta, flowing from left to right, for the case of Eq. ( |2^ ) 
in Fig. |H (C). We note that the presence of the (non- 
perturbative) instanton-induced interaction results in 
an integration over the momentum exchanged by the 
quarks through the instanton, similar to the presence of 
perturbative gluon exchanges leading to loop integrals, 
Fig. |2| (B). This analogy still holds when multi- instanton 
effects are taken into account: each instanton gives an 
additional "loop" integral over the instanton induced mo- 
mentum transfer. Furthermore, we read off the "Feyn- 
man rules" : each free quark propagator contributes an 
exponential and each zero-mode a modified Bessel func- 
tion. Combined, the W2W three-point function, Eq (||), 
in the SIA reads 



n _ ^(frcc) 
Lt4 — yjr^ 



G 



(C) 



G 



(D) 



(28) 



Time reversal invariance imposes that the contributions 
from Figs. ^ (C) and (D) are equal, as can be seen from 
Eqs. ( |26| ) and (|27| ) by substituting t — > — t and q — » q 
(recall the absolute value in the exponential of the free 
W2W propagator, Eq. @). Finally, the W2W two- 
point function, Eq. (|l0|), in the SIA, is given by 



G z 



dp I dziii(p)- 



i*2 



-|q|/2 (f(* 4 -2t)+£(24)) 



J-oo 4m* 
(l + |q|£(* 4 -2t)/2) (l + |q|e(z 4 )/2) 



£(z 4 - 2tf 



and analogously 



G = G (frcc) + G z 



(29) 



(30) 



The pion formfactor is now readily obtained by means 
of Eq. ([)]). At the distances we are considering (t = 
0.7 fm), the free contribution to the correlation functions, 

G^ Tee ^ and G^ hee \ are at most 30% corrections (at |q| = 
5 GeV) to the zero-mode three- and two-point functions. 
This can be seen by evaluating e.g., the zero- mode two- 
point function for large |q|: 



G zm — 



dp 



n{p) p 4 
16 m* 2 , 



I 2tt 



| q |3/2 e -|q|«t) 



(31) 



5 When one considers many instanton effects, the use of the singu- 
lar gauge is obligatory, as one desires the topological charge to 
be localized. In singular gauge, however, the expression of the 
zero-mode W2W quark propagator is more complicated. 



by a saddle point analysis of the z 4 integration. In com- 
parison to the free contribution, Eq. (|22|), this leads to a 
30% correction at |q| = 5 GeV, for standard values of the 
average instanton density n = 1 fm~ 4 and the quark effec- 
tive mass to* = 70MeV. We consider such a correction 
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to be the edge of the window we work in 6 . Therefore, 
we may approximate the formfactor by 7 



U-4 



(t a- 



(C) + (D)^ q. 



-*,-f) 



G( frec )(2i,q/2) + G, m (2i,q/2) 



G 



4 2 m 



G 2m (2i, q/2) 



(32) 



Consequently, the effective parameter m* and the average 
instanton density n, which are multi-instanton induced 
model parameters, are irrelevant for the pion formfactor 
at intermediate momentum transfer. We stress that, in 
the approach presented here, we do not need to evalu- 
ate the pion wave function, since we obtain directly the 
relevant Green functions from the SIA. In other words, 
the pion wave function is implicitly included in our cal- 
culation through the large time extraction of the pion 
pole. 



III. NUMERICAL RESULTS 

Before comparing our results to the experimental data, 
we check that the ratio of the zero-mode three-point to 
two-point function, Eq. (|32|), does not depend on the 
distance t, which ensures that the pion pole contribution 
has been successfully isolated. For this purpose, we show 



n — 1 — 1 — 1 — r 



■ ■ ■ t = 0.1 fm 
-- t = 0.6fm 

— t = 0.7fm 

— Mon. Fit 
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FIG. 3: The ratio of three- to two-point function, Eq. (p2[), 
evaluated in the SIA for different values of t and for simplicity 
a delta function distribution at p = 1/3 fm. At t > 0.5 fm our 
results are t independent. 



6 Note, however, that similar corrections appear in the numerator 
and the denominator, and thus the real accuracy of the formfac- 
tor can actually be better. 

7 We note that the free correlators do not have a pion pole. They 
are the leading order 0(a®) continuum contributions. The rel- 
evance of these terms at higher momentum transfer demands 
larger times t. As discussed, however, then the SIA breaks down. 



this ratio for several distances t in Fig. |[ We find that 
the pion pole contribution dominates for t > 0.5 fm, in 
agreement with |l4|, [27) . 

Our central result for the pion formfactor is plotted in 
Fig. Q in comparison to the JL AB data and the monopole 
formfactor fitted to the low Q 2 data. We observe that 
the SIA prediction agrees with the recent JLAB mea- 
surements. In addition, it is completely consistent with 
the monopole form in the kinematic region accessible to 
JLAB. These results complement the analysis of Blotz 
and Shuryak [fl9|| , where the same agreement was found 
at small momentum transfer. Upcoming measurements 
at JLAB will be able to test the single instanton predic- 
tion as the microscopic mechanism at intermediate mo- 
mentum transfer. 

We find that the instanton contribution remains well 
above the pQCD scale set by Eq. (|l]) throughout the 
experimental region under investigation at JLAB, Q 2 < 
6 GeV 2 . 

As the momentum is increased, one requires larger 
times in order to isolate the pion pole from the higher 
excitations. In sum-rule approaches, contributions from 
the continuum are obtained from the free and perturba- 
tive gluon exchange correlation functions, e.g., Fig. ^ (B). 
Perturbative contributions will never develop a pion pole. 
Indeed, we observed that the free correlators, G^ 00 "* and 
G (frco) , become non negligible for |q| > 4 - 5 GeV. The 
presence of such continuum contributions destroys the 
time-independence of the ratio of three- to two-point 
functions around t ~ 0.6 fm, shown in Fig. ^, and hence 
the method can no longer be used. 

In principle, one could extend the range of validity of 
our approach by increasing the time. However we recall 
that, when t becomes larger, multi-instanton effects are 
significant and the SIA breaks down. We conclude that 
our approach is not applicable to the study of the high 
momentum transfer region, Q 2 > 20 GeV 2 . 

Let us now discuss, how one can get the transition 
to the pQCD limit from instantons. For large momen- 
tum transfer, where counting rules and factorization are 
justified, the evaluation of the formfactor requires non- 
perturbative information, encoded in the pion wave func- 




(E) 



(F) 



FIG. 4: (E): Diagram contributing to the formfactor at 
asymptotic momentum transfer in pQCD. The blobs de- 
note the pion wave functions <f>. (F) A leading order multi- 
instanton contribution to the W2W three-point correlator re- 
sponsible for the transition toward the pQCD result. 
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tion c6, in addition to pQCD diagrams. This is shown 
diagrammatically in Fig. |l| (E). In the instanton liquid 
model, the pion is bound because the quark and anti- 
quark feel a strong attraction from instantons. As we 
argued, at large momentum transfers, one has to incor- 
porate multi-instanton effects in order to isolate the pion 
pole and the pion wave function. For example, we give 
a leading order multi-instanton diagram relevant in the 
asymptotic limit in Fig. [I| (F). 

Since the pion formfactor may experimentally be 
rather accurately measured, it is instructive to ask 
whether such data may shed some light on the instan- 
ton size distribution. In Fig. [jj, we have plotted the 
results of our theoretical predictions for Q 2 F 7T (Q 2 ) ob- 
tained for different cases of n(p). We contrast the sim- 
plest size distribution Jig , n{p) ~ S(p — 1/3 fm), to the 
results obtained from a lattice QCD parameterization. 
We notice that the presence of tails in the size distri- 
bution introduces only small corrections to the formfac- 
tor. Therefore, we conclude that the simplest distribu- 
tion n(p) ~ 5(p — 1/3 fm) indeed captures the relevant 
features for the pion formfactor at intermediate momen- 
tum transfer. We observe that our result becomes closer 
to the perturbative limit, if the average instanton size is 
larger or possibly if there is an asymmetric tail toward 
larger-sized instantons in the distribution. 

Next, we compare the results from distributions which 
have the same small-size limit, but cut off at different in- 
stanton sizes, p — 0.37 fm and 0.47 fm. We observe that, 
throughout the entire kinematic region we have consid- 
ered, our predictions are quite sensitive to the average 
instanton size. This implies that the asymptotic region, 
where the pion formfactor probes the small-size 't Hooft 
behavior, is not reached within the window, where our 
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FIG. 5: The dependence of the pion formfactor Q 2 F 7r (Q 2 ) on 
the instanton size distribution. The SIA (solid) curve repre- 
sents a small-size 't Hooft distribution with a Fermi distribu- 
tion cutoff and lattice QCD parameters. The SIA II (dotted) 
curve has a different mean instanton size p = 0.47 fm (same 
width) and the SIA III (dashed) curve is obtained with the 
simplest delta distribution n(p) = n5(p — 1/3 fm). 



approach is justified. 

Summarizing, we have shown that the presence of en- 
hanced instanton-induced forces constitutes a large and 
clearly dominant contribution to the charged pion form- 
factor. The latter is so large that the perturbative regime 
will not be reached experimentally. This nicely contrasts 
the situation for the 77*71"° neutral pion transition form- 
factor, where the asymptotic pQCD regime is reached 
much earlier, at Q 2 ~ 2 GeV 2 |§. This striking dif- 
ference is explained by instanton arguments as well, be- 
cause the two corresponding three-point functions have 
a different chiral structure. Physically, this is the same 
reason why the vector and axial channels have a rather 
strong "Zweig" rule, forbidding flavor mixing, while for 
the pseudoscalars such a mixing is very strong. 



IV. CONCLUSIONS 

We have presented the first calculation of the instanton 
contribution to the charged pion electro-magnetic form- 
factor in the momentum range Q 2 > lGeV 2 . Techni- 
cally, it is based on the ratio of three- to two-point cor- 
relators, in which the enhancement comes from the zero- 
mode quark W2W propagators. Because it is based on a 
representation in momentum space, the framework devel- 
oped presents several quite direct analogies with ordinary 
perturbative calculus. 

Our calculation shows that, when the non-perturbative 
instanton-induced forces are taken into account, the 
charged pion formfactor remains much larger than the 
asymptotic perturbative prediction, throughout the en- 
tire experimentally accessible region of JLAB. This result 
is in contrast to other model calculations, which do not 
include such a force @, § §, [u], [ll], ||, || . Clearly, the 
upcoming JLAB measurements will provide a unique op- 
portunity to test the role played by the instanton-induced 
't Hooft interaction in hadronic physics. 

Moreover, we found that the SIA prediction coincides 
with the monopole form up to Q 2 ~ 10 GeV 2 . We em- 
phasize that this does not represent a test of the validity 
of our model, since the monopole fit is a phenomenolog- 
ical parametrization exclusively of the low energy data. 
Nevertheless, such an agreement is interesting, because 
it suggests that the widely used vector dominance model 
for the pion is in fact a dynamic property of the QCD 
vacuum, and that its validity extends into the interme- 
diate momentum regime. 

We finally studied whether our results depend on the 
details of the instanton size distribution. We found that 
the pion formfactor is determined by the average instan- 
ton size up to rather large Q 2 and is insensitive to the 
details of the size distributions. 

In addition, an application of the same method to the 
analysis of the proton electric formfactor has been car- 
ried out and leads to a dipole formfactor at intermediate 
momentum transfers fl32|. 
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